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Abstract 

The ground and thermal states of a quadratic hamiltonian representing the interaction of bosonic 
modes or particles are always Gaussian states. We investigate the entanglement properties of these 
states for the case where the interactions are represented by harmonic forces acting along the 
edges of symmetric graphs, i.e. 1, 2, and 3 dimensional rectangular lattices, mean held clusters and 
platonic solids. We determine the Entanglement of Formation (EoF) as a function of the interaction 
strength, calculate the maximum EoF in each case and compare these values with the bounds found 
in £Q which are valid for any quadratic hamiltonian. 
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1 Introduction 



Suppose we have a many body quantum system in its ground state. We ask how much quantum 
correlation or entanglement exists between two of the particles? How this entanglement depends on 
the strength of interaction between the particles? Does it extend far beyond the nearest neighbor 
particles? If we raise the temperature, hence mixing the ground state with higher level states, at 
what temperature the entanglement ceases to exist? Some of these questions have been investigated 
in recent years for spin systems j^J El IU El 13 El ■ Quite recently another class of systems composed of 
bosonic modes of systems with continuous degrees of freedom have come under intensive investigations 
El ini El El QH El El El El- There are a number of reasons for this interest. First, bosonic 
modes are the appropriate subsystems the entanglement of which must be calculated when dealing 
with systems of identical boson particles ^1 EOl ESj . Second, states of continuous systems are 
widely encountered in many branches of physics in which entanglement plays a role, i.e. namely in 
quantum optical setups, in atomic ensembles interacting with electromagnetic fields (§], in the motion 
of ions in ion traps and in low excitations of bosonic field theories. Third, for a class of such states, 
namely Gaussian states, analytical measures of entanglement have been defined and calculated in 
closed form [15] . Finally, a large class of interesting many body systems of this type, when written in 
terms of suitable coordinates, are in fact free systems. This latter reason makes a fuller investigation of 
the above questions in such systems much easier than in spin systems. For example it has been shown 
|12j that in the ground state of a ring of particles coupled with each other by harmonic forces, entan- 
glement exists only between adjacent particles and is exactly zero if the particles are non-adjacent. It 
has also been possible to study the development of entanglement when such systems evolve in time 
|13j . Neither of these questions are easy to investigate in spin systems. 

Furthermore quite interestingly the notion of entanglement frustration has recently been introduced in 
PQ . The motivation for this notion comes from the observation that quantum entanglement like many 
other local properties of a system becomes strongly restricted when one imposes global symmetry 
requirements on the state of a system of N particles. For example the Entanglement of Formation 
(EoF) ^21 of two modes in a Gaussian state is generally unbounded, it can assume any value between 
zero and infinity. However if these modes are part of a three mode system having a permutation 
symmetry, then the amount of their EoF becomes finite. Therefore one may ask how the EoF of two 
modes in a system, may depend on the total number of modes, the symmetry of the state and the 
dimension of the lattice which represents the interaction of these modes. These questions have been 
studied in , with the aim of determining the maximum value of EoF which can exist between two 
modes given that the whole state has a general symmetry. However the Hamiltonian of is a specific 
hamiltonian having no tunable parameter (i.e. frequencies), in fact it has been so defined that its 
ground state allows the maximum value of EoF among all the possible quadratic Hamiltonians with 
the same symmetry. 

In this paper we want to study the effect of interaction strength and symmetry in determining the 
entanglement of bosonic modes in a class of states pertaining to the ground state of symmetric graphs. 
We take a system of N particles or bosonic modes coupled quadratically with each other so that the 
Hamiltonian allows a global symmetry. We then determine the entanglement of two modes when the 
system is in its ground state. The entanglement depends on the number of particles, the strength of 
interaction and the symmetry of the Hamiltonian. We will see that entanglement always increases 
(almost rapidly) with the strength of interaction and always saturates to a finite limiting value but 
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does not exceed the bounds found in pQ. We also show that in some cases symmetry has a constructive 
effect on entanglement. For example in some symmetric graphs like a cube or an octahedron, there are 
entanglement between next-nearest neighbors which is known to be absent in less symmetric graphs, 
i.e. in rings |EsJ. In the most simple case, where the graph consists of two vertices connected by an 
edge, we also calculate the thermal entanglement as a function of temperature and interaction strength 
and determine the threshold temperature beyond which entanglement vanishes. 

The structure of this paper is as follows: In section |2 we review briefly the definition of Gaussian 
states and the Entanglement of Formation (EoF) for symmetric two-mode Gaussian states. In section 
131 we introduce the quadratic Hamiltonian and calculate the EoF of two Gaussian modes when the 
whole system is in the ground state of this quadratic hamiltonian. This hamiltonian simply describes a 
mass-spring system of the kind considered in |T2] and is different from the one introduced in pQ , which 
maximizes the EoF. In section 0] we present our results for various symmetric graphs. We end the 
paper with a conclusion. A remark on notation. Throughout the paper we use E^ n f for the maximum 
EoF which is obtained at infinite frequencies (interaction strength of the springs) in our hamiltonian 
and E max for the bound discovered in pQ. 



2 Gaussian States 



In this section we discuss the rudimentary material on Gaussian states |23| that we need in the sequel. 
Reference [23] can be consulted for a rather detailed review on the subject of Gaussian states. 
Let R := (xi,X2, ■ ■ ■ xjy,pi,p2, • • -Pn) be N conjugate operators characterizing ./V modes and subject 
to the canonical commutation relations 



where a 



In 



[Rk,Ri] = i(?ki, 

is the 2n dimensional symplectic matrix and I n denotes the n dimensional 



unit matrix. 

A quantum state p is called Gaussian if its characteristic function defined as 

C(0 :=tr(e-^ fc<Tfci %), 
is a Gaussian function of the £ variables, namely if 

C(f) := e ^~ ikTkli \ 

where we have assumed that linear terms have been removed by suitable unitary transformations. The 
matrix T, called the covariance matrix of the state, encodes all the correlations in the form 



kl ■' 



(R k R l + R l R k )-2{R k ){R l ). 



By symplectic transformations the covariance matrix of a two mode symmetric Gaussian state (one 
which is invariant under the interchange of the two modes) can always be put into the standard form 



/ n k 



V 



\ 

k p n J 



(1) 
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where k x > > k p and k x > \k p \. For such a state the Entanglement of Formation (EoF) has been 
obtained in |15) . It is given by 

EoF(p) := C+ log 2 C + - C_ log 2 C_ , (2) 

in which 



C± = (3) 



and 



A := min(l, 6 := J(n — k x )(n + kp)). (4) 

Thus a state is entangled only if 8 < 1. In general -Eoi* 1 is unbounded and can assume any value 
between for 5 = 1 to infinity for 5 — ► 0. 

3 The ground and thermal entanglement of a system of particles 
coupled by harmonic forces 

Consider now a system of particles with canonical variables R = (x±,X2, ■ ■ ■ XN,pi,p%, ■ ■ -Pn) T subject 
to the following Hamiltonian 

1 N 1 N 1 

k=l k=l k,l 

where the prime over the sum indicates that only specific particles are coupled with each other. This 
Hamiltonian can be written in matrix form 

*-iW?jW (a) 



where V introduces the quadratic matrix of the potential. The state of such a system at temperature 
T = -s is given by p = 6 z , where Z = tre~^ H is the partition function. It is not difficult to show 
that this state is Gaussian with covariance matrix 



r — © T p , 

T x = W^coth^p 
BW 

T p = T^coth^-, (7) 

where W := V?. At zero temperature when only the ground state is populated the covariance 
matrix tends to 

r(r = 0) = W' 1 © W. (8) 

Once the covariance matrix of the iV modes is obtained as above, the covariance matrix pertaining to 
any two modes, say the modes 1 and 2 is given by the sub- matrix obtained from V by deleting the rows 
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and columns corresponding to all the other modes. Furthermore if this Gaussian state is symmetric, 
i.e. if it is invariant under the interchange of the two modes, then its covariance matrix will have the 
following (not-yet standard) form 
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n p 
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where 



n x 


= (r x ) u 


= (r.; 


22 


n p : 




= (T P ) 


22 


m x 


:= (r x )n 


= (r x 


)21 


rn p 


■■= (r P )i2 


= (r P 


)21 



By the canonical transformation 



1 

-Pl,2, 

a 



%1,2 >■ OiX lt 2, Pi, 2 ~ 

where a = the covariance matrix Q takes the standard form 

\ 
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( n k x 
k% in 

V 



n k 

k 
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n ) 



where 



(9) 



(10) 
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(12) 



n 
k x 

kn 



m. 



rrir 



n 3 



(13) 



From this standard form and (JU one can now determine the EoF of the two modes. 
Inserting (fl3|) in Q we find: 



which in view of (fTU|) leads to 



(n - k x )(n + k p ) 



I n 



n. 



{yjn x n p - m x J- E -){^n x n p + m v A — ) 



n 



n,, 



(n x - m x )(n p + m p ), 



5 = y/tfTxhi - (r,)i 2 ) ((r p )n + (r p ) 12 ). 



(14) 



(15) 
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The matrix elements of T x and T p are obtained by diagonalization of V. At zero temperature (i.e. for 
the ground state) we will have 
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^((W-^n - (W- 1 )^) (W u + W 12 ). 



(16) 



In the following section we will apply these results to evaluate the EoF of states for modes pertaining 
to various symmetric lattices. 

4 Entanglement in symmetric graphs 

Let G be a symmetric graph, having TV" vertices. The adjacency matrix of the graph is denoted by A, 
where by definition A^ = 1 if the the vertices i and j are linked by an edge and Aij = otherwise. 
In a concrete situation we can think of the particles of unit mass located at the vertices of the graph 
experiencing a global harmonic potential and coupled through springs of strength k = to 2 with each 
other. If the number of nearest neighbors of a vertex is denoted by z, then the potential matrix V 
takes the following form 



In the following we will consider various symmetric graphs and obtain in each case the EoF of the 
two modes corresponding to various pair of vertices. In almost all cases calculations show that en- 
tanglement exist only between nearest neighbor sites. The only exceptions are the octahedron and 
the cube where the next-nearest and next-next-nearest neighbors are also entangled with each other. 
In the simplest case, where the lattice consists of only two vertices connected by an edge, we also 
determine the threshold temperature above which entanglement disappears. This case may model the 
entanglement between the vibrational modes of two atoms in a molecule or two ions in an ion trap 
which is destroyed by thermal fluctuations. 

4.1 The simplest two- vertex graph 

The simplest symmetric graph consists of two vertices connected by an edge, figure (^-a). The potential 
for this graph is V = \{x\ + x\ + co 2 (xi — X2) 2 ). It represents the harmonic interaction of two atoms 
with frequency u each trapped inside a harmonic well of unit frequency. The potential matrix is 




(17) 




(18) 



with eigenvalues and eigenvectors 




(19) 



This leads to the following values for the relevant matrix elements: 
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Figure 1: The simplest symmetric graphs. The graph in (b) is only symmetric with respect to the 
interchange of the vertices 1 and 3. The EoF in graph (a) is unbounded, that of graph (c) is bounded 
but much higher than that of vertices 1 and 3 in graph (b). 



The equation 5 = 1 determines the threshold temperature. Inserting this value of delta in © and 
(J5J) will determine the EoF as a function of the frequency u and the temperature, shown in fig- 
ure ©. The maximum entanglement obtains at zero temperature and is unbounded, it behaves as 
E = —0.56 + log 2 (w) + 0{uj~ 2 ) for large frequencies. The threshold temperature increases almost 
linearly with frequency. 

4.2 The simplest three-vertex symmetric graphs 

We now consider the graph shown in figure (^-b). This is one of the two symmetric graphs with three 
vertices. The other graph which is a triangle will be treated as a special case of iV-vertices mean field 
clusters (for N = 3) . Note that this graph is symmetric only under the interchange of vertices 1 and 
3. We will determine the EoF of these two modes. Our motivation for studying this graph is to show 
that entanglement can exist between next-nearest neighbors in open chains with an odd number of 
vertices, a property which we have seen also for higher than 3 vertices in our calculations. Here the 
potential is 




(20) 



Inserting these values in 1)15(1 we find the value of 5, 




(21) 
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T 



Figure 2: Thermal entanglement of the two modes (particles) (in units of ebits) associated to the 
graph in figure (JDa), as a function of temperature and frequency. E is dimensionless and we are 
working in units in which co and T are also dimensionless. 



V = o ( X l + X 2 + ^3) + -7r(( x l ~ x 2) + {X2 ~ X 3 



(22) 



with the potential matrix 



V 



( 1 + co 2 

V 



The eigenvalues and eigenvectors are 



l07 



1+uJ 2 , 



—UJ 

1 + 2uj 2 

-UJ 2 



|ei) 

l e 2> 
l e 3) 





-uj 2 

1+UJ 2 



^(1,1, If, 

1 T 

7i (1 '- 2 ' 1) 



(23) 



(24) 



In this case we are only interested in the ground state entanglement. So only the matrix elements of 
W and W~ l need be calculated. One finds after straightforward calculations along the lines in the 
previous subsection the following value for 8, 



I2 + V1 + 3uj 2 
3V1 + oj 2 



(25) 



The value of this parameter is always less than 1, and hence there is entanglement between the 
next-nearest neighbors at all frequencies. The interesting point is however that this entanglement is 
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bounded. It obtains at very large frequencies where 5 — ► 3 * , leading to Ei n f = 13.62 in units of 
10" 2 ebits. 

4.3 Mean field clusters 

Consider a mean field cluster of N vertices in which every vertex is connected to N — 1 other vertices. 
The potential matrix for this graph is given by 

V = (l + Nlo 2 )I - cu 2 E, (26) 

where E is the matrix all of whose entries are equal to 1, Eu = 1 V i, and j. Using the property 
E 2 = NE, we find 



1 - Vl + Nu 2 



W = Vl + Nu 2 I+ — — E 

Using these matrices and equation (|16|) one finds the value of 5 at zero temperature 



S N = J 2 + (iSr - 2)VIT ^ (28) 

Inserting this into equations © and © gives the entanglement as a function of N and frequency. 
Let us first consider the special case of 3 vertices, whose graph is a triangle, shown in figure (JDc) and 
compare it with result for the three- vertex graph shown in figure (|flb). 

For N = 3, we have 5 3 = yj 2 ^^* - The Eof saturates at E inf = 40.08 (in units of 10~ 2 ebits) for 
infinite frequencies. 

Returning to (|28[) we can find the EoF in mean field clusters as a function of frequency for fixed number 
of vertices. The maximum EoF always occurs for infinitely large frequencies where 5 tends to x ~ J 



N ■ 

Figure ) shows the EoF as a function of w for some mean field clusters with N = 4,6, 8, 12 and 
20 vertices. These are the number of vertices of regular polyhedra, namely tetrahedron, octahedron, 
cube, dodecahedron and isocahedron respectively, which will be treated in the nest subsection. 



Finally figure© shows how the maximum EoF decreases with increasing the number of vertices in 

■2io g -- 

2W 



mean field clusters. For large values of N the EoF vanishes as E ~ 1,72 ^ ogN 



4. 4 Regular Polyhedra 

In this section we consider the regular polyhedra, the tetrahedron (N = 4), the cube [N = 8), the 
octahedron (N = 6), the dodecahedron (N = 20) and the isocahedron [N = 12). The tetrahedron is 
a mean field cluster with N = 4, so we can use the results of the mean field cluster for this graph. For 
the other cases we follow the same procedure as above, i.e. determine the adjacency matrix which in 
turn determines the potential matrix from Q17|). The square root W and the inverse square root W^ 1 
of V are calculated by diagonalization of V. We can then calculate the necessary elements of W and 
W -1 in order to determine the EoF of pair of vertices in terms of the strength of the interaction u. 
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w 

Figure 3: The EoF (in units of 10 -2 ebits) of two vertices for mean field graphs with number of vertices 
from top to bottom N = 4, 6, 8, 12 and 20. 

It turns out that the EoF is always an increasing function of uj and saturates for infinite frequencies. 
The maximum EoF attains its value for the nearest neighbors. Depending on the polyhedron, the EoF 
may or may not exist between next nearest neighbors. Furthermore compared with mean field clusters 
the EoF is always higher than that of mean field clusters, the difference is larger for higher frequencies. 
Figure © shows the EoF of the regular polyhedra in terms of the strength of interaction. Figure (jHJ) 
compares the EoF of the octahedron with that of a mean field cluster with 6 vertices. Finally table 
(1) compares the maximum EoF obtained in regular polyhedra with those of the mean field clusters 
of the same size (the same number of vertices) and the bounds found in pQ. 
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Polyhedron 


N 


Mean field Cluster 




E 

-'-'max 


Tetrahedron 


4 


19.74 


19.74 


19.74 


Cube 


8 


8.30 


9.80 


19.74 


Octahedron 


6 


4.68 


9.74 


10.75 


Dodecahedron 


20 


2.15 


7.00 


11.12 


Isocahedron 


12 


0.83 


4.51 


5.37 



Table 1-The maximum entanglement (Ei n f) in regular polyhedra and mean field graphs of the same 
size, compared with the bounds found in [Q. N is the number of vertices of the polyhedron and the 

mean field graph. 
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Figure 4: The maximum EoF (in units of 10 2 ebits) for mean field clusters as a function of the 
number of vertices. 
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0.24 






Octahedron 


9.74 


2.58 








Dodecahedron 


7.00 














Isocahedron 


4.51 















Table 2-The EoF between different vertices in regular polyhedra. The superscript indicates the 
minimum distance between the vertices in terms of the number of edges. The — sign indicates that 
such distances do not occur in the corresponding regular polyhedron. 

4.5 Rings and rectangular lattices 

For the d-dimensional cubic lattice, of N d vertices, the potential matrix is 



a 

v = (i + idu?)i - cu^in + n' 1 ), 



(29) 



k=i 



where is the cyclic unit shift operator in the k— direction. The eigenvalues and eigenvectors are as 
follows: 

^ 2 = l+4 W 2 Esin 2 ^, Us (n):=-^=e^. (30) 



k=i 



10 




w 



Figure 5: The EoF of adjacent vertices (in units 10 2 ebits) in regular polyhedra. The curves corre- 
spond from top to bottom to Tetrahedron, Cube, Octahedron, Dodecahedron and Isocahedron. 

From these values one can easily obtain the relevant matrix elements of W and W~ l . For example 
for the one dimensional lattice we will have 

N-i 

W ±l k ,i= £(l + 4u; 2 sin 2 ^)~e^r^. (31) 

s=0 

As always the EoF increases with the strength of interaction. It turns out that EoF exists only between 
nearest neighbor sites. Figure (j7J) compares the EoF of rings with two very different size as a function 
of frequency. The interesting point is that, the difference for the number of vertices only shows itself 
when the frequency is high enough. 

Figure (JHJ) compares the EoF of d— dimensional lattices for d = 1,2,3 when N — > oo. For the infinite 
lattices, the sums in Q31|) turn into integrals and we will have for two adjacent vertices say 1 and 2 of 
one dimensional lattice, 

Wtx = - r dx(l + 4uj 2 sin 2 x)^ 

7T Jo 

W^2 = ~ r dx< y l + 4uj2 sin2 cos 2x ' ( 32 ) 

7T Jo 

and for the two dimensional lattice 

W^ 1 = I dx dy[l + 4o; 2 (sin 2 a; + sin 2 y)] ± 2 
tt Jo Jo^ 

W^ 1 = — / dx I dy[l + 4a; 2 (sin 2 x + sin 2 y)] ± ^cos23;. (33) 

7T Jo Jo 

Table (3) compares the saturation values of EoF achieved in these cases with the maximum achievable 
EoF's derived in pQ. 
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Figure 6: The EoF of adjacent vertices (in units 10 2 ebits) as a function of w for the regular 
octahedron (line) and a six- vertex mean field cluster (dots). 



d 




E 


1 


16.34 


30 


2 


3.54 


6.31 


3 


1.50 


2.62 



Table 3-The maximum EoF between adjacent vertices in d-dimensional rectangular lattices 

compared with the bounds found in pQ. 

5 Discussion 

We have investigated the entanglement properties of the ground state of a quadratic Hamiltonian 
governing the harmonic interaction of bosonic modes or particles located on the vertices of symmetric 
graphs i.e. 1,2, and 3 dimensional rectangular lattices, mean field clusters and platonic solids. The 
entanglement has been calculated as a function of the interaction strength. In each case the EoF 
is an increasing function of the interaction strength and it saturates to a finite value dictated by 
the type of the graph. Our results confirm those of in that the maximum values that we obtain 
are always less than the bounds found in that paper which correspond to a particular quadratic 
hamiltonian proved to allow for maximum achievable EoF between any two modes in the class of all 
quadratic hamiltonians. Some peripheral results may be interesting. First, in some symmetric graphs 
( the cube, the octahedron and open chains with an odd number of vertices) there are entanglement 
between non-adjacent vertices, a property which is extraordinary for entanglement, since for rings it 
has been shown that entanglement can not extend beyond the nearest neighbors. Second, in the 
case of two atoms vibrating in a molecule or ions vibrating in an ion trap, corresponding to the simple 
graph (Jfla) the threshold temperature above which entanglement is destroyed, has been determined. 
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EoF 



Figure 7: The EoF of two adjacent sites for rings with different number of vertices almost coincide for 
low frequencies, N = 11 (points) N = 201 (line). 
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